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Single-particle spectral properties near the Mott transition in the one-dimensional Hubbard model 
are investigated by using the dynamical density-matrix renormalization group method and the Bethe 
ansatz. The pseudogap, hole-pocket behavior, spectral-weight transfer, and upper Hubbard band arc 
explained in terms of spinons, holons, antiholons, and doublons. The Mott transition is characterized 
by the emergence of a gapless mode whose dispersion relation extends up to the order of hopping t 
(spin exchange J) in the weak (strong) interaction regime caused by infinitesimal doping. 
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PACS numbers: 71.30.+h, 71.10.Fd, 74.72.Kf, 79.60.-i 

High-T c cuprate superconductors are obtained by dop- 
ing Mott insulators that have gapless spin and gapped 
charge excitations. In the small-doping regime, anoma- 
lous features, such as pseudogaps, Fermi arcs, hole pock- 
ets, and kinks in dispersion relations, have been observed 
[l|, l2|. It is widely believed that understanding such 
anomalous electronic properties near the Mott transition 
will be critical for achieving high-T c superconductivity. 

The Mott transition can be viewed as charge localiza- 
tion caused by strong Coulomb repulsions; note that the 
nature of this transition contrasts with that of conven- 
tional metal-to-band-insulator transitions where a band 
of single-particle states is fully occupied regardless of the 
interactions. Further, because of the strong correlations, 
it is generally difficult to obtain reliable information and 
intuitive understanding on the Mott transition. How- 
ever, although some detailed properties might depend on 
the lattice structure, it can be expected that essential 
features of Mott transitions are generally true for finite 
dimensions and that the generic features can be deduced 
by using a one-dimensional (ID) Hubbard model. In this 
model, we can interpret excitations relevant to the Mott 
transition using exact solutions without bias. Also, we 
can investigate the properties of a single metallic phase 
all the way to the Mott transition point without causing 
any instability to phase separation, superconducting, or 
magnetic orders, which might otherwise have occurred in 
higher dimensions. 

In this Letter, we investigate the spectral properties 
near the Mott transition in a ID Hubbard chain using un- 
biased numerical techniques and exact solutions, focusing 
attention on the pseudogap, hole-pocket behavior, upper 
Hubbard band (UHB), and spectral- weight transfer from 
the UHB to the lower Hubbard band (LHB). We also dis- 
cuss the nature of the Mott transition through compar- 
isons with the Fermi liquid picture and recent numerical 
results on the two-dimensional (2D) Hubbard model [3[. 

Model and method.— We consider the ID repulsive 
Hubbard model defined by the following Hamiltonian: 



H 



■ 7 / \ C i+lc 



H.c 



•"£< 



n-t n n ~ /-* 



£-< 



where Ci a and m a are the annihilation and number oper- 
ators of an electron at site i with spin a, respectively, and 
Ui = riif + mi- The hopping integral t and on-site repul- 
sion U are positive. The number of sites, electrons, and 
down spins is denoted by L, N, and M, respectively. The 
doping concentration is defined as 5 = 1 — N/L. Also, 
we define the single-particle spectral function as follows: 



A(k,u >0) 
A(k,u <0) 
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Here, |GS) and \l) denote the ground state with energy 
£"gs and the excited state with Ei, respectively, while 
c, creates an electron with momentum k and spin a. 

I calculated A(k,ui) by using the dynamical density- 
matrix renormalization group (DDMRG) method [4{ un- 
der the open boundary condition with the number of 
density-matrix eigenstates m = 120. The data with 
Lorentzian broadening with half width at half maximum 
r\ = 0. 16t were deconvolved to those of Gaussian broaden- 
ing with standard deviation a = O.lt. The spectral func- 
tions at k = ttj/(L + 1) for j = 1 ~ L were extrapolated 
to k — and w. Figure [1] (a) shows the typical behavior 
of A(k,u>). Since the difference between the results for 
L = 60 and 80 is small in this scale, this Letter shows the 
results for L — 60 unless otherwise mentioned. Noting 
that A(k, uS) = A(—k,cj) = A(k + 2tt,ui), we consider the 
properties of < k < tt without loss of generality. 

To identify the dominant modes in A(k,oj), we use 
the Bethe ansatz, where the wave functions, energies, 
and momenta of eigenstates are expressed in terms of 
{kj} and {A Q } that satisfy the following Bethe equa- 
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= 1 - M. Here I j 
are integers or half-odd integers. We impose 
the (anti)periodic boundary condition, when M in the 
ground state is odd (even). Because eigenstates are ob- 
tained through the Bethe equations once {Ij} and {J a } 
are given, the {Ij} and {J a } distributions characterize 
eigenstates. In the ground state, {Ij} and { J a } are con- 
secutively distributed around zero. Excited states are ob- 
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FIG. 1: (a) A(k,oo) for L=60 (red solid line) and 80 (blue 
dotted line) at fc~0.8?r for U/t-8 when 5=0.1. (b) A(u) 
in the LHB (red lines) for (b-1) U/t=8, 5~0.133; (b-2) 
U/t=&, <5~0.033; and (b-3) U/t=4, (5~0.033. In (b-1), the 
blue dashed-dotted line denotes the contribution of dominant 
modes, while the green region shows the contribution of con- 
tinua except the dominant modes. In (b-2, 3), the light blue 
regions denote the contribution from uj>0. Arrows indicate 
e(k=n) of Eq. (JTJ. (c) Dispersion relations of the upper edge 
of the sh* continuum (red dotted line) and the h* mode with 
s(kp) (brown solid line) for U/t—8 at <5~0.133. The green 
dotted line and the dot at fe=0.57r show those of S— 5*0. (d) 
A(k)=f du;A(k y uj) for U/t=8 at <5~0.133. Blue squares with 
the light blue region denote the contribution from uj<0, which 
is the momentum distribution function in the ground state 
n(k). Green diamonds denote the contribution from oj>0; 
red crosses and pink open circles denote that of the LHB for 
w>0 and that of the UHB. The vertical yellow line indicates 
k—kF- Brown solid circles denote the total weight at each 
k, which satisfies the sum rule within numerical accuracy, (e) 
Doping dependence of spectral weight A. Solid lines show A of 
the LHB for uj>0 at U/t—i (red circles) and 8 (blue squares). 
Red and blue triangles denote A for uj<0 at U/t—A and 8. 
The green dashed line and orange dashed-dotted line indicate 
A— 5 and (1— <S)/2, respectively. Open symbols denote data 
for L=80. (f) A(k,uj) for J7/t=4 at <5~0.033. 



tained by creating holes (particles) inside (outside) the 
consecutive distributions. The holes in the {J a } distri- 
butions are called spinons, while those in the {Ij} distri- 
butions are called holons. The particles created outside 
the consecutive {Ij} distribution are called antiholons. 
Hereafter, we use the following shorthand notations for 
the spinon, holon, and antiholon: s, h, and h*, respec- 
tively. Noting that the momentum of the eigenstate is 



given by K = 2-7r(V ■ Ij + J2 a Ja)/L, we define the mo- 
menta of s, h, and h* as q s = 2irJ s /L, q^ = 2irIh/L, 
and qh* = 2irIh*/L, respectively, where J s , 1%, and Ih* 
are the positions of s, h, and h* in the distributions of 
{J a } and {Ij}. The momentum ranges are \q s \ < kp, 
\qh\ < 2k p, and 2k p < \qh*\ < tt with the Fermi momen- 
tum k F = tt(1 - S)/2. 

Dominant modes for lo < 0.— Before discussing the 
Mott transition, we review the dominant modes for lo < 
at 5 = 0.4 @, 0j- The dominant mode for k < k F = 0.3n 
near w < in Fig. [2](a-l) is identified as the s mode with- 
out h and h* [blue dashed-dotted line for k < kp in Fig. 
[2jc-l)]. The mode slightly below it is identified as the h 
mode with s having q s — kp and h* having gv = 2kp. 
Hereafter, we denote it as the h mode with s(kp) and 
h*(2kp). The mode for k < 3kp = 0.97T is identified as 
the h mode with s(—kp) and h*(—2kp), which is called 
the shadow band [8(- The lower edge of the shh* contin- 
uum is indicated by the lowest pink dotted line in Fig. 
[2^c-l). For lo > 0, sh* excitations are dominant [7|, |9(. 

Pseudogap.— As shown by the red line in Fig. QJb- 
1), the momentum-integrated spectral weight [.A (a;) = 
/- l§^(' c ' a -')] shows reductions near lo ~ 0. This pseu- 
dogap behavior originates from the following three prop- 
erties: (i) low-energy property as a Tomonaga-Luttinger 
liquid lfj, [llj , (ii) the band-edge singularity of the ID 



dominant modes, and (iii) continua spread above and be- 
low us — 0. In a Tomonaga-Luttinger liquid, A(w) oc \oo\ e 
for lo — > with the same exponent as that of the mo- 
mentum distribution function n(k) — 1/2 oc \k — kp\ e for 
k — > kp, where < 6 < 1/8 depending on the values of 
U/t and 6 [Fig. Hfd)] [jj|. Because 9 -> 1/8 as 8 -> 0, 
the pseudogap behavior will be significant near the Mott 
transition. Also, the contribution of dominant modes ex- 
tracted by Gaussian fitting of the peaks in A(k, to) at each 
k [blue dashed-dotted line in Fig. [TJb-1)] shows peaks 
near the bottom of the s mode for lo < and the top of 
the h* mode with s(kp) [brown solid line near k ~ ir/2 
in Fig. [2fc-2)] for lo > 0; this implies that the pseudogap 
behavior can be explained as a dip between the peaks 
near the band edges of the ID dominant modes. More- 
over, the contribution from continua except the dominant 
modes [green region in Fig. [TJb-1)] shows reductions near 
lo ~ 0. As in Fig. [2Jb-2), the continua above and below 
lo = shrink to the gapless points at k — kp and 2ir—3kp 
as lo — > 0. The continua for lo > are mainly due to the 
2-A-string solutions where two A a 's are A ± iL T /(4t) with 
real A [12| for k < kp, sh* for k ~ 27r — 3kp, and hh* 
excitations with s(kp) and h*(±2kp) for k ~ 2ir — 3kp 
and k F . For lo < 0, the continua mainly come from the 
sh excitations with h*(±2kp). 

Hole-pocket behavior.— The mode connecting the two 
gapless points at k = kF and 27r — 3&f for w > is 
identified as the h* mode with s(kp). On the basis that 
it has charge character and that the momentum region 
between these gapless points shrinks as Afc = 2ir5, this 



5=0.4 



5^0.133 



5=0.033 




^ 


(e-2)" ■ ' 


(e-3)- ■ ' ■ 



0.5 

kin, 



0.5 

kin 



0.5 

kin 



FIG. 2: (a),(b),(d) A(k,u) for U/t=8 at <5~0.4, 0.133, 0.033 
(from the left) for (a) overall views, (b) the LHB, and (d) the 
UHB. (c),(e) Dispersion relations obtained using the Bethe 
ansatz in L=240, corresponding to (b),(d). Solid lines except 
those at uj=0 show holon modes and antiholon modes. Blue 
dashed-dotted lines denote spinon modes. Red (pink) dotted 
lines show upper (lower) edges of continua of hh* with s(fcp) 
and h*(2kp), sh* , hh* with s(kp), and shh* (from above) 
in (c), and shd in (e). Purple dashed double-dotted lines in 
(c) denote spinon modes in 2-A-string solutions. Open circles 
indicate peaks of dominant modes of (b),(d). The inset in (c- 
3) is the close-up near the gapless points. The pink [yellow] 
region in (e) denotes the hd continuum with s( — Uf) [ s (^f)]- 
Spinon [holon] modes in (e) have h{±2kp) and d (±(-7T — 2/c_f)) 
with q d qh<0 [s(±k F ) and d(|g d |~7r<5/2, q d qh<0)]. 



region can be regarded as a hole pocket. It is expected 
that the h* mode will become robust as U/t increases 
because of the reduction in double occupancies. In fact, 
the h* mode has been found for U/t — > oo [8| and in the 
ID t-J model O]. 

Spectral-weight transfer.— Although the spectral 
weight A transferred to the LHB for uo > equals the 
amount of doping (A = S) at t = [14j, Fig. OJe) shows 
that A > S for t ^ 0. Such behaviors have been discussed 



in the literature H 0. [ For w < 0, ^4 = (1 - c$)/2 
as usual.] The question here is which mode carries the 
transferred spectral weight. As in Figs. HJb-3) and (c-3), 
a considerable weight is carried by the mode of the up- 
per edge of the sh* continuum. In the 5 — > limit, where 
\qh'\ — > 7r, this mode reduces to the s mode with h*(±n); 
the dispersion relation e(k) is obtained as follows: 



e =-57- /_„. dq cos 2 q sech 
k =l + ITJZ dx J~ir d( l sech 
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by using dressed energies and momenta for S — > 16lll7|. 
This gapless mode naturally leads to the 2-spinon con- 
tinuum at half-fillin g exp ressed as E(k) = e(fei) + 6(^2) 
with k — k\ + k 2 — 7T [16|, [l7| : The two gapless modes car- 
rying 5=1/2 with gapless points at k = ±kp cause the 
gapless 5=1 spin excitations with the gapless point at 
|fc| = 2kF- It should be noted that the cosine dispersion 
relation in the U/t — )• limit as well as the flat dispersion 
relation in the U/t — > 00 limit [18J are reproduced, since 
e(k) ^ -2tcosk for U <t and e(k) 
J 

plies that the spectral weights transferred by infinitesimal 
doping spread up to the energy of 0(t) for small U/t and 
O(J) for large U/t, as in Figs. QIb-2)-(c), (f), and^b- 
3) . This behavior contrasts with band-insulator-to- metal 
transitions, where the transferred spectral weight re- 
mains within the chemical potential shift A/j, of 0(8 2 ^ D ) 
for S — > in D dimensions. 

Upper Hubbard band.— To identify the dominant modes 
in the UHB, we consider the solutions with a k-A string, 
i.e., a pair of complex kj's with a nonzero imaginary part, 



Tp cos k with 
4t 2 /U for U > t (H| for \k ± n\ < tt/2. This im- 



which represents a pair of electrons 12| . Since such so- 
lutions have energies of 0(U) in the large U/t regime 
[121 Il6j . it is natural to expect that they are relevant 
for the UHB [9|, [l9(. In fact, it has been found that 
solutions with a string length greater than 1 have con- 
siderable spectral wei ght s in the high-energy regime for 
quasi- ID spin models 2fJ,|21| and spinless-fermion chains 
22J, |23j . The k-A string is characterized by a half-odd in- 
teger or an integer I\. We call the quasiparticle (QP) for 
the k-A string a doublon and denote it as d; the momen- 
tum is defined as q^ = 2wIi/L in |g<j| < ir—2kp = irS [12J | . 
We define the sign of q^ such that k = —q s — qh+Qh* +1d- 
At half-filling, A(k,u>) of the UHB is symmetric with 
that of the LHB because of the particle-hole symmetry. 
In the small-J regime, A(k,u>) remains almost symmet- 
ric, as in Figs. [TJf) and[2ja-3). In the 5 — > limit, where 
\qd\ — > 0, the dominant modes in the UHB are character- 
ized by s and h in the fc-A-string solutions and have es- 
sentially the same character as that in the LHB. Namely, 
the low-energy mode for k > kp ~ 7r/2 in Fig. [UJd-3) 
is identified as the s mode with h(±2kp) and d{qd — 0), 
as in Fig. EJe-3). The mode slightly above it and the 
low-energy mode for k < kp are both identified as the h 
mode with s(kp) and d(qd — 0). The high-energy mode 



is identified as the h mode with s(—kp) and d(qd ~ 0), 
which corresponds to the shadow band. For larger dop- 
ing, the dominant modes with charge character in Figs. 
[DJd-1,2) remain within the hd continua with s(±kp) [yel- 
low and pink regions in Figs. [He- 1,2)] and appear to have 
d{\q d \ ~ ir5/2, q d q h < 0) [solid lines in Figs, ^e-1,2)}. 

Discussion.— To obtain an intuitive understanding of 
Mott transitions, two scenarios have been considered. 
One is that the effective carrier density (n*) vanishes, as 
in the case of metal-to-band-insulator transitions. The 
other is the effective-mass (m*) divergence in the Fermi 
liquid theory. The present results show that spectral 
weights are mainly transferred to the h* mode with s(fcp) 
and the mode of the upper edge of the sh* continuum 
[Figs. QXb-1)— (c)]. The former behaves like n* — > 0, 
because the hole pocket shrinks as 5 — > 0. The latter 
may seem to behave as m* — > oo, because its spectral 
weight decreases [Fig. [He)]. However, in a Fermi liquid, 
the renormalized bandwidth near the chemical potential 
shrinks as e(fc) oc 1/m* for m* — > oo, which is not the 
case for the present results wherein e(k) in Eq. (JTJ re- 
mains 0(t) or O(J) [Fig. Die)]. Thus, a simple classifica- 
tion of n* — » or m* — > oo is not accurate for the Mott 
transition of the Hubbard chain; the Mott transition is 
rather characterized as a loss of charge character from 
the mode having both spin and charge characters, while 
the dispersion relation of the spin part remains gapless 
and dispersing. 

The nature of the Mott transition in the Hubbard chain 
may be more clearly understood through comparisons 
with the results of a recent cellular dynamical mean- 
field study on the Mott transition of the 2D Hubbard 
model in Ref. [3j , where the following features have been 
suggested, (i) The phase transition from a large Fermi 
surface to hole pockets occurs at nonzero S. (ii) The 
transferred spectral weight behaves as A = S + (ni^ni\) 
after tiny doping, (hi) A gap fully opens between the 
band forming the hole pockets and that of ingap states 
for u > below the UHB, leading to the zero surface of 
the Green function [24] . (iv) The pseudogap can be at- 
tributed to doublon (doubly occupied site)-holon (empty 
site) binding [15| relaxed due to doping. 

In contrast, for the Hubbard chain, (i) there is no phase 
transition when 5^0: The gapless point at k = 2n — ikp 
persists even in the large doping regime, where the gap- 
less point goes beyond the zone boundary (fc = 7r) and is 
located at k = 3kp, as in Fig. EJb-1). (ii) The transferred 
spectral weight does not follow A = 6 + (n^ni^) and 
gradually reaches zero as 8 — > 0, as in Fig. QJe). (iii) The 
cosinelikc mode and the antiholon mode are both gapless 
at k — kp: There is no gap between them at k = kp. 
The gapless cosinelike mode is related to the gapless spin 
excitations, (iv) Double occupancy in the LHB does not 
behave as a QP. The QP for a pair of electrons is defined 
in the fc-A-string solutions, which are relevant only for 
the UHB: The pseudogap behavior in the LHB is not re- 



lated to the QP. Understanding and explaining the above 
differences can be an interesting topic for future study. 

Summary.— The single-particle spectral properties 
near the Mott transition were investigated by using the 
Bethe ansatz and DDMRG method in the ID Hubbard 
model. Characteristic spectral features near the Mott 
transition, such as the pseudogap, hole-pocket behav- 
ior, spectral-weight transfer, and upper Hubbard band, 
were explained in a unified manner in terms of spinons, 
holons, antiholons, and doublons. A remarkable feature 
is the emergence of the gapless mode extending up to 
0(t) [O(J)] for small [large] U/t by infinitesimal doping. 
This mode is related to the spin excitations at half-filling 
and cannot be interpreted by either the rigid-band pic- 
ture or the Fermi liquid theory. The present results have 
similarities to the anomalous spectral features observed 



in high-T c cuprates 
Hubbard model [3j, [l 




and some aspects of the 2D 
24j . I expect that generic features 
of Mott transitions in finite dimensions can be deduced 
from the present results. 
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